#include <omd.h>

class MyMDClass:public MDSystem {

    void CreateSystem() {
        AddAtom(new FCC("111",30,20,10, "platinum"))
           ->Create()
           ->SetTemperature(100.0)
           ->SetName("Crystal");
    }

    void CreateGadget() {
        SetIntegrator(new MDIntegrator);
        AddForce(new TForceEAM("platinum"));
        AddModify(new VerletList);
        AddDetector(new SysMonitor("md.out"));  
        AddDetector(new ThermoDetector(0.02));
    }

    void SystemSetting() {
        PBoundary=NONPERIODIC;
        MaxTime=0.2;
        SetOutputDirectory("output");
    }

    void BeforeRun() {
        DumpAtoms("init.dat");
        PrintInfo("info.out");
    }

};

int main(int argc, char* argv[]) {
    MyMDClass TheSim;
    TheSim.SetArgument(argc,argv);
    TheSim.Run();
    return 0;
}
